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A GLOBAL EMPIRICAL MODEL OF 

THERMOSPHERIC COMPOSITION BASED ON 

OGO-6 MASS SPECTROMETER MEASUREMENTS 

A. E. Hedin 
H. G. Mayr 
C. A. Reber 
G. R. Carignan* 

N. W. Spencer 

A global empirical model for magnetically quiet conditions has been devel- 
oped to describe longitudinally averaged OGO-6 N 2 , O, and He densities (in the 
altitude range 400 to 600km) by means of an expansion into spherical harmonics. 
The accuracy of the analytic description is of the order of the experimental er- 
ror, thus providing an excellent representation of the satellite observations. 

Some of the results are summarized. (1) The annual asymmetric component in 
N 2 density corresponds to an increase of about 400°K between winter and sum- 
mer poles, which is about a factor of three greater than in the Jacchia model. 
However, the corresponding oxygen component at satellite altitudes reflects an 
increase between winter and summer comparable to that in the Jacchia model; 
this implies that at 120km atomic oxygen increases by more than a factor of two 
between summer and winter poles. (2) The amplitude of the symmeti'ic semi- 
annual component in N 2 is larger at the poles than the equator. At the poles the 
N 2 variation corresponds to a temperature decrease of about 70°K from mid-April 
to mid-July. (3) At low latitudes the daily temperature maximum occurs close to 
16 hours local time, nearly in agreement with radar backscatter observations. 

♦Space Physics Research Laboratory, University of Michigan, Ann Arbor. Michigan. 
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A GLOBAL EMPIRICAL MODEL OF 
THERMOSPHERIC COMPOSITION BASED ON 
OGO-6 MASS SPECTROMETER MEASUREMENTS 

Until recently, the bulk of the data available for determining the global be- 
haviour of neutral density in the thermosphere has consisted of total density 
measurements obtained from satellite drag. Variations in composition can not 
be uniquely determined from these data and the atmospheric models' ’ 7 based 
on satellite drag have, until recently, generally assumed diffusive equilibrium 
with fixed boundary densities at 120km. 

This paper reports on some aspects of the global behaviour of N 2 , O, and 
He as determined from data obtained with the neutral mass spectrometer carried 
aboard the OGO-6 satellite. This satellite was launched on June 5, 1969 into an 
82° inclination orbit with a 398km perigee and 1100km apogee. The data cover 
the time period June 27, 1969 to June 22, 1970 during which time data were ob- 
tained in the altitude range between perigee and 600km (800km for He). In order 

to study the behaviour of the thermosphere independent of longitude and universal 
(3 4) 

time effects' ’ 7 the globe was divided into 5° latitude bands (with a further sub- 
division near the poles into south and north bound passes) and data falling within 
each band were averaged over a 24 hour (universal time) period. Data from 
days which had an Ap index greater than 5 were excluded in order to minimize 
magnetic activity effects. ^ 

Two different schemes have been used to describe the data. First, the 
longitudinal averages have been fit with a formula based on a spherical harmonic 
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expansion (in geographic latitude - local time coordinates) in which the composi- 
tion is described in terms of density variations at 450km with an altitude depend- 
ence assuming isothermal conditions at any given time. The expansion includes 
local time independent, diurnal, and semidiurnal terms and the coefficients have 
been written as the sums of terms representing the annual, semiannual, solar 
flux, and magnetic activity variations. The temperature used in the calculation 
was determined from the N ? data in a manner to be described shortly, but it is in 
any case not critical to determining 450 km densities as this altitude is near the 
average altitude of the measurements. A total of 33 parameters was determined 
by a least squares fit to the data for each constituent. The average error of the 
fit was 20% for N 2 (about twice the experimental error) and 10% for He and O 
(about equal to the experimental error). 

A second approach used to describe the measured densities was analogous 
to the method employed in the Jacchia J65 static diffusion model. With the 
density of N 2 fixed at 120km, the exospheric temperature necessary to calculate 
the measured densities was determined using analytic density profiles suggested 

tr* \ 

by Walker . ' ’ The temperature derivation was based on the N 2 density, rather 
than O or He, because this species is theoretically closest to being in diffusive 

(7) 

equilibrium. However, the average temperature needed to fit the measure- 
ments is not arbitrary, as the measurements are made at different altitudes, 
and there is an average temperature which will fit the data best. Details of the 
global variation of temperature cannot, however, be obtained with much signifi- 
cance from the scale height information. The procedure followed was to fix the 
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120km N 2 density, fit the measured densities using a spherical harmonic for- 
mula for exospheric temperature variation analogous to that used for the 450km 
density variations, and allow the temperature gradient at 120km to be an addi- 
tional free parameter so that the average exospheric temperature and density 
would be consistent with the measurements. By using the model temperatures 
determined from N 2 , the average densities and variations at 120km needed to 
reproduce the measurements of O and He were determined. It must be kept 
clearly in mind that although this model represents the measured data quite well, 
the densities at 120km are only inferred on the basis of the diffusive equilibrium 
assumption. However, the variations in composition at 120km are a convenient 
way of illustrating departures from fixed boundary diffusive equilibrium models 
and to the extent that diffusive equilibrium holds, the O to N 2 and He to N 2 ratios 
inferred at 120km are relatively independent of errors in the derived exospheric 
temperature. 

. The global average densities at 450km (average F 107 = 150, Ap = 3. 5), . 

determined as described above, are 2.52 x 10 6 cm -3 for N 2 , 6.42 x 10 7 cm' 3 
for O, and 3. 55 x 10 6 cm -3 for He. Assuming an N 2 density at 120km of 3 x 
10 11 cm -3 and a temperature of 355°K the inferred densities at 120 km are 
8.81 x 10 10 cm* 3 for O, and 3.20 x 10 7 cm -3 for He. The average exospheric 
temperature was 1116°K. The S value in the temperature profile^ was 0.021. 

In discussing particular density variations, the results will be given in terms 
of departures from the above averages. 
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DIURNAL AND SEMIDIURNAL COMPONENTS 


The variation of density with local time at the equator for equinox conditions 
is shown in Figure 1. At 450km the N 2 density has a maximum at 15.6 hours, 
the O density at 14. 9 hours and the He density at 10.2 hours. The phase of the 
O density maximum is not significantly later than that in the J65 model. As 
compared to the J65 model, the later temperature phase and larger maximum 

/ Q \ 

to minimum ratio is in agreement with incoherent scatter observations. The 

amplitudes of the N 2 and O diurnal density variations are in agreement with the 

J65 model and the small (statistically insignificant) inferred O density variation 

at 120km indicates that the measured 450km N 2 and O densities are reasonably 

consistent with a fixed boundary diffusive equilibrium model. The morning 

maximum for He densities and the systematic change in phase from N 2 to O to 

He, however, indicate that the thermosphere may not be in complete diffusive 

equilibrium and the measured trend is consistent with predictions of dynamic 

(7) 

effects due to winds . 

ASYMMETRICAL ANNUAL COMPONENTS 

The seasonal variations in density lead to the latitudinal density profiles 
seen in Figure 2 for solstice conditions. The measured oxygen density variation 
at 450km is comparable to that in the J65 model, but the N 2 density variation is 
about three times larger. Thus, contrary to the situation for the diurnal varia- 
tions, the relative variation of N 2 and O densities shows a striking departure 
from a fixed boundary diffusive equilibrium model and the inferred 120km oxygen 
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densities show a winter density enhancement of 2.4 times over summer values. 
The possibility of a winter oxygen enhancement has been inferred from iono- 
spheric data^ 9, and a similar trend has been recently deduced from a com- 
bination of drag and mass spectrometer data. The greater exospheric tem- 
perature variation, from winter to summer pole, deduced from the N 2 densities 
is in the right direction to resolve a recent' ' objection to explaining the F- 
region seasonal anomaly by changes in the oxygen to nitrogen ratio. 

(13, 14) 

The measured helium densities clearly reflect the winter helium bulge, 

however, the inferred helium density variation at 120km is about double that 

(2) 

given in a recent Jacchia model. v ' The latitude profile in the summer hemi- 
sphere is not as well defined by the OGO-6 measurements as in the winter hemi- 
sphere because of a lack of sensitivity, but the indicated gradient is consistent 

( 15) 

with measurements from Explorer 32. 

The latitudinal-seasonal variations in both helium and oxygen are consistent 
with a circulation pattern in the thermosphere which has air rising near the 
summer pole and transporting air rich in helium and oxygen to the winter 
pole.' 16 ' 17 ' 18 ' 19) 

SYMMETRICAL SEMIANNUAL COMPONENTS 

The latitudinal profile of the amplitude of the semiannual variation is illus- 
trated in Figure 3. The amplitudes are quite small and are believed to be ac- 
curate to no more than 5 or 10%. Nevertheless, the fact that the phase of the 
variation was independently determined for each gas to be mid- April (or 
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mid-October) leads to some confidence that the semiannual effect has indeed been 
detected. Furthermore, the latitudinal variations of the amplitudes show a sys- 
tematic trend from nitrogen to oxygen to helium. The oxygen amplitude shows 
only a slight variation with latitude, consistent with J65. The increase of N 2 
densities toward the poles (and presumably temperature) combined with the 

relative depletion of oxygen and helium is qualitatively consistent with dynamical 
(19) 

models.' ' The importance of the composition is also in line with the recently 
revised Jacchia model^^ which represents the semiannual effect in terms of 
density rather than temperature variations. 

SUMMARY 

The analysis of composition measurements made with the neutral mass 
spectrometer aboard the OGO-6 satellite leads to the following conclusions. The 
measured atomic oxygen densities are generally in good agreement with those 
deduced from drag. The molecular nitrogen densities in the annual and semi- 
annual variations depart significantly from those predicted by drag models and 
suggest similar departures for exospheric temperatures. The helium densities 
generally tend to vary in an inverse manner to the nitrogen densities. These 
composition changes are consistent with dynamical processes associated with 
the global circulation in the thermosphere. 
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FIGURE CAPTIONS 


1. Diurnal and semidiurnal components of density and temperature variation 
from the OGO-6 empirical density model for the equator at equinox. 

2. Asymmetrical annual component of density and temperature variation from 
the OGO-6 empirical density model for solstice conditions. 

3. Symmetrical semiannual component of density and temperature variation 
from the OGO-6 empirical density model. 
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